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ABSTRACT 

We examine the Kepler light curves of VI 504 Cyg and V344 Lyr, encompassing ~460 d at 1 min cadence. 
During this span each system exhibited ~40 outbursts, including four superoutbursts. We find that, in both 
systems, the normal outbursts lying between two superoutbursts increase in duration by a factor ~ 1.2— 1.7, 
and then reset to a small value after the following superoutburst. In V344 Lyr the trend of quiescent intervals 
between normal outbursts is to increase to a local maximum about half way through the supercycle — the in- 
terval from one superoutburst to the next — and then to decrease back to a small value by the time of the next 
superoutburst, whereas for V1504 Cyg the quiescent intervals are relatively constant during the supercycle. 

Both of these trends are inconsistent with the Osaki’s thermal-tidal model, which robustly predicts a secular 
increase in the quiescent intervals between normal outbursts during a supercycle. Also, most of the normal out- 
bursts have an asymmetric, fast-rise/slower-decline shape, which would be consistent with outbursts triggered 
at large radii. The exponential rate of decay of the plateau phase of the superoutbursts is 8 d mag -1 for VI 504 
Cyg and 12 d mag -4 for V344 Lyr. This time scale gives a direct measure of the viscous time scale in the 
outer accretion disk given the expectation that the entire disk is in the hot, viscous state during superoutburst. 

The resulting constraint on the Shakura-Sunyaev parameter, a hot — 0.1, is consistent with the value inferred 
from the fast dwarf nova decays. By looking at the slow decay rate for superoutbursts, which occur in systems 
below the period gap, in combination with the slow decay rate in one long outburst above the period gap (in U 
Gem), we infer a steep dependence of the decay rate on orbital period for long outbursts. We argue that this 
relation implies a steep dependence of a co id on orbital period, which may be consistent with recent findings of 
Patterson, and is consistent with tidal torquing as being the dominant angular momentum transport mechanism 
in quiescent disks in interacting binary systems. 

Subject headings: cataclysmic variables: dwarf novae 


1. INTRODUCTION 

Cataclysmic variables (CVs) are semi-detached interacting 
binaries in which a Roche lobe filling K or M secondary trans- 
fers matter to a white dwarf (WD). CVs evolve to shorter or- 
bital periods and show a “gap” between Ph = 2 and 3 (where 
Ph = -Porbitai/1 hr) during which time the secondary star 
loses contact with its Roche lobe and mass transfer ceases. 
Thus the binary becomes fully detached. At Ph — 2 the 
secondary comes back into contact with its Roche lobe and 
mass transfer resumes. For Ph < 2 angular momentum loss 
from the binary is thought to be due solely to gravitational 
radiation. The CV subclass of dwarf novae (DNe) are char- 
acterized by their semi-periodic outbursts. SU UMa stars are 
DNe lying below the period gap that exhibit short, normal 
outbursts (NOs) and superoutbursts (SOs). We refer the time 
from one SO to the next as the supercycle. SOs show su- 
perhumps which are modulations in the light curve at periods 
slightly exceeding the orbital period. There are two further 
subdivisions within the SU UMa grouping: (i) the VW Hyi 
stars at long orbital periods, near Ph — 2, for which the decay 
rate is fairly constant during a SO, and (ii) the WZ Sge stars at 
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short orbital periods, a little greater than Ph = 1, which have 
less frequent, larger amplitude SOs, for which the decay rate 
decreases during a SO. DNe outbursts are thought to be due to 
a limit cycle accretion disk instability (Lasota 2001) in which 
material is accumulated in quiescence and then dumped onto 
the WD during outburst. During short outbursts in longer pe- 
riod DNe, a few percent of the stored mass is accreted, and 
during long outbursts a significant fraction ~0.2 of the stored 
mass is accreted. For the SU UMa stars, a SO is thought to 
accrete > 0.7 — 0.8 of the stored mass. Although the accretion 
disk is never in steady state during the limit cycle, it is close to 
steady state during SO, with the local rate of accretion inside 
the disk \M(r)\ decreasing linearly from a maximum value at 
the disk inner edge to zero at the outer edge. The accretion 
disk modeling has traditionally been done within the Shakura 
& Sunyaev (1973, hereafter SS) formalism, using two values 
for the a viscosity parameter, ahot for gas in the hot, ionized 
disk, and a co id for gas in the quiescent disk. 

There are two bright SU UMa stars in the Kepler field ex- 
hibiting a variety of temporal behavior that make them worthy 
of a detailed statistical study of their outbursts, VI 504 Cyg 
(Kepler ID 7446357; P h = 1.67) and V344 Lyr ( Kepler ID 
7659570; Ph = 2.1). These are members of the VW Hyi sub- 
division. To date the two light curves have amassed 459.8 d at 
1 —min cadence. Excluding gaps and bad data points, the light 
curves contain 637948 and 637915 entries, respectively. Pre- 
vious studies of the Kepler data on SU UMa stars have found 
quiescent superhumps in V344 Lyr (Still et al. 2010), pre- 
sented numerical models of the long term light curve of V344 
Lyr (Cannizzo et al. 2010; hereafter CIO), and studied super- 
humps, both positive and negative, in the long term V344 Lyr 



light curve (Wood et al. 201 ! ). 

Statistical studies of DNe have been useful in delineating 
the long-term properties of outbursts, e.g.., the outburst du- 
ration, recurrence times, and quiescent intervals, and placing 
physical constraints on models of outbursts (e.g., Campbell 
& Shapley 1940, Sterne, Campbell & Shapley 1940, Bath 
& van Paradijs 1983, van Paradijs 1983, Szkody & Mattei 
1984, Cannizzo & Mattei 1992, 1998, Ak, et al. 2002, Simon 
2004). Several interesting studies of the SU UMa stars have 
been carried out. For instance, van Paradijs (1983) compared 
the variation in outburst duration 4 for long and short out- 
bursts above the 2 — 3 hr CV period gap with those below, and 
found that long outbursts appear to be more significant in con- 
trast to short outbursts as one progresses to shorter orbital pe- 
riod systems only because the durations of the short outbursts 
scale with orbital period, but those of the long outbursts are 
relatively constant. Therefore the ratio of outburst durations 
f&(SO)/f&(NO) is large. This finding was amplified by Ak et 
al. (2002). In this work we seek to extend previous studies 
by analyzing the high fidelity Kepler light curves of two SU 
UMa stars to investigate the properties of the outbursts. In 
Section 2 we examine the outburst properties of the NOs and 
SOs, in Section 3 we discuss the results, in particular the scal- 
ing of the SO decay rate with orbital period, and in Section 4 
we summarize our results. 

2. OUTBURST PROPERTIES 
2.1. Normal Outbursts 

The relatively small number of total outbursts observed in 
the two systems to date does not yet allow one to amass de- 
tailed frequency histogram distributions, as was done for in- 
stance by Cannizzo & Mattei (1992) for the ~700 outbursts 
observed up to that time in SS Cyg. However, the variations of 
the outburst duration tb and intervening quiescent interval t,, 
within a supercycle contain information which can potentially 
guide theoretical models. For instance, in the thermal-tidal 
model for the supercycle in the SU UMa systems (Ichikawa & 
Osaki 1992, 1994; Ichikawa, Hirose, & Osaki 1993), one sees 
a monotonically increasing sequence of NO quiescent inter- 
vals leading up to a SO. Now with the exquisite Kepler data, 
this can be directly tested. The NOs of the SU UMa stars have 
been poorly characterized prior to Kepler due to their intrinsic 
faintness. 

Figures 1 and 2 present the currently available long term 
light curves for VI 504 Cyg and V344 Lyr. To obtain the 
apparent Kepler magnitude from the original data, electron 
flux per cadence, a correspondence 7 of Kp — 12 to 10 7 e - 
cadence -1 was adopted. Each light curve encompasses four 
SOs; VI 504 Cyg had 36 NOs while V344 Lyr had 33. Two 
systematic effects are worth noting: (i) There are some gaps 
in the data due to monthly data downloads and quarterly 90° 
spacecraft rolls, (ii) In addition, the long-term absolute flux 
calibration for Kepler can be problematic for several reasons 
(Haas et al. 2010. Koch et al. 2010): (1) Target masks do not 
contain all the flux from the target. The pixel mask is chosen 
to optimize signal-to-noise for each target. This is done for 
the mid-time of the quarter. As Kepler orbits the sun. how- 
ever. targets drift across the mask as a result of differential 
velocity aberration. Varying amounts of target flux fall out- 
side of the pixel aperture over the quarter, so the target flux 
tends to vary smoothly over the quarter due to this effect. (2) 

' See http://keplergo.arc.nasa.gov/CaIibrationSN.shtml. 


After each Kepler rotation, the target falls on a different CCD 
chip with a different PSF and the target is centered differently, 
e.g., perhaps on a pixel boundary rather than a pixel center. 
The 4 arcsec pixels undersample the PSF. This means that the 
pixel mask changes shape and size each quarter, with differ- 
ent amount of light losses. (3) With different shaped masks, 
the number and amount of contaminating sources in the mask 
will differ across the quarter gap. 

Problematic long term flux calibration is evident in the ele- 
vated flux level seen in V1504 Cyg during 100 < t{ d) < 200 
(i.e, Q2 of observations). In addition there can also be a sec- 
ular trending of the flux calibration within a quarter, which 
may account for the apparent ~50 d e— folding decrease of 
the quiescent level following the first two SOs in V344 Lyr, 
both of which occurred just after a roll maneuver. In view 
of these considerations, we limit our attention in this study to 
timescales rather than amplitudes. 

In a standard light curve plot of magnitude versus time, a 

straight line decay means that the flux / = /oe -t// ' r , which 
we would therefore refer to as an exponential decay. A de- 
cay which deviates from exponential could be concave up- 
ward (slower-than-exponential), concave downward (faster- 
than-exponential), or some more complicated shape. The 
decays of the NOs are faster-than-exponential, which obvi- 
ates theoretical studies that attempt to introduce a weak radial 
power law dependency into the SS a parameter in order to get 
precisely exponential decays (e.g., Ichikawa, Hirose, & Osaki 
1993; Cannizzo 1994). We may quantify the deviation from 
exponentiality in the NO decays by dividing the decay into 
bins, each spanning a small range A Kp = 0.15 — 0.20, and 
calculating the locally defined decay rate td(d mag -1 ) within 
each bin. Dividing the decay into a total of 10 bins for each 
star, we take A Kp = 0.15 for V1504 Cyg and A Kp = 0.20 
for V344 Lyr. For demonstration purposes, the binning dis- 
tribution for the decay of the first NO in V34 Lyr is shown 
in Figure 3. Within each bin we compute a frequency his- 
togram distribution of td values for all NO decays. The re- 
sults are shown in Figures 4 and 5. The small numbers in each 
panel indicate the magnitude range in Kp measured from the 
local maximum for each NO. Thus, the upper panels show 
the decay rates computed near maximum light, and as one 
progresses to lower panels one is sampling the fainter por- 
tions. For both systems one sees a shifting in the centroid 
of the frequency histogram distribution to lower td values as 
the decay proceeds, equivalent to saying that the decays are 
concave downward. In both systems, the parts of the decays 
fainter than about 1 mag below NO maximum asymptote to a 
centroid median of ~0.6— 0.7 d mag -1 . Thus the concave- 
downward aspect of the NO decays is due in large part to 
the rounded maxima; by the time the NO has faded to within 
A Kp ss 1 mag of quiescence, the rest of the decay is closely 
exponential. 

Bailey (1975) discovered that the decay rate for NOs is cor- 
related with orbital period (see his Figure 1 for eight well- 
studied systems), and derived the empirical relation td — 0.38 
d mag -1 Ph. Smak (1984) was the first to show that this re- 
lation can be used to infer «hot — 0.1. Subsequent work 
confirmed and expanded on this finding (Cannizzo 2001a). 
Warner ( 1 995b) looked at a larger sample than Bailey and de- 
duced the fit td = 0.53 d mag -i P^ 84 . For the two systems of 
interest in this study, for which P h = 1.67 (V1504 Cyg) and 
2. 1 (V344 Lyr). respectively, the predictions for the rates of 
decay would be 0.64 and 0.80 d mag -1 from the original Bai- 
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FiG. 1 . — Kepler light curve for VI 504 Cyg at 1 min cadence, spanning ~460 d. The light curve covers four superoutbursts and 36 normal outbursts. The 
vertical red lines indicate the local maxima for the normal outbursts in which coverage permits a reliable determination. 
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Fig. 3. — The distribution of decay bins for the first normal outburst in 
V344 Lyr. 

ley relation, and 0.82 and 0.99 d mag -1 from the Warner fit. 
The presently determined values ~0.7 d mag -1 in V1504 Cyg 
and ~0.6 d mag -1 in V344 Lyr for the faintest bins are closer 
to Bailey’s original expression. If one averaged over the entire 
decay, the values would be larger, closer to Warner’s expres- 
sion. The precision of the current determinations are limited 
by the small number of total NOs in both systems. With time, 
better statistics will allow better determinations, as has been 
done for instance with SS Cyg (Cannizzo & Mattei 1998, see 
their Figure 6). 

Figure 6 shows the outburst durations f;, for the two sys- 
tems. Based on the magnitude range for outbursts in the 
two systems, we use a cut of Kp — 16 in V1504 Cyg and 
Kp = 17.9 in V344 Lyr to define outburst onset and termina- 
tion. There is an obvious trend of increasing outburst width 
tb between consecutive SOs: during a supercycle, increases 
by 20—70%. Figure 7 shows the asymmetry parameter f r , 
defined to be the ratio of the rise time to the total outburst 
duration, for the individual NOs depicted in Figures 1 and 
2. Figure 8 shows the quiescent intervals t q between succes- 
sive outbursts. For most of the time, the quiescent intervals 
in VI 504 Cyg are relatively constant. For the last two super- 
cycles there is an increase in t q following a SO. In contrast, 
for V344 Lyr the t q values increase to a local maximum value 
about half way between successive SOs, and then decrease 
back to roughly their starting values by the time of the next 
SO. 

At roughly the time of the third SO in VI 504 Cyg, strong 
negative superhumps appear (Wood et al. 201 1, in prepara- 
tion) and persist through the end of the light curve. The pres- 
ence of negative superhumps indicates that the mass is not 
accreting primarily at the rim of the disk. Although there is 
not an obvious change in f r at this time, it is perhaps note- 
worthy that t q increases. This may indicate a strong drop in 
the mass transfer rate from the secondary star, which would 
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Fig. 4. — A histogram of the fast decay rates for the 36 normal outbursts in 
VI 504 Cyg with well-defined maxima, as shown in Figure 1. The numbers 
in each panel give the magnitude range in A Kp down from maximum cor- 
responding to the interval within which the local decay rate 4 is measured, 
thus brighter portions of the decay are in the upper panels. 


TABLE 1 

Summary of Parameters 


Parameter 

VI 504 Cyg 

V344 Lyr 

Comment 

NO: 4(d) 

1. 1-3.0 

2.4— 4.9 

NO duration 

SO: 4(d) 

11.9 

6.8 

SO duration 

NO: f r 

0.18-0.44 

0.22-0.37 

NO rise/duration 

i«(d) 

5.3-17.6 

2.4-21.4 

Quiescent duration 

NO: t,j(d mag l ) 

0.7 

0.6 

NO decay rate 

SO: 4(d mag" 1 ) 

7.9 

12 

SO decay rate 


increase the recurrence time for outbursts. In the magnetic 
systems such as AM Her, which lack accretion disks, one has 
direct evidence for the instantaneous mass transfer rate from 
the secondary stars in the long term light curves, and one does 
often see intervals during which the mass transfer rate reduces 
to a low level (e.g., Hessman, Giinsicke, & Mattei 2000). 

2.2. Superoutbursts 

Figure 9 shows a composite of the four SOs in V 1 504 Cyg, 
positioned in time so that t — 4 d corresponds to the crossing 
of Kp — 16 during SO onset. The rate of decay during SO is 
~8 d mag'" 1 . Figure 10 shows a composite of the four SOs in 
V344 Lyr, positioned in time so that t ■ 3 d corresponds to 
the crossing of Kp = 18 during SO onset. The rate of decay 
during SO is ~12 d mag" 1 . Within each system the SOs are 
quite similar in overall appearance and duration. 

3. DISCUSSION 
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Fig. 5. — A histogram of the fast decay rates for the 28 normal outbursts in 
V344 Cyg with well-defined maxima and well-sampled decays, as shown in 
Figure 2. Conventions are the same as in Figure 4. 

A summary of the parameters associated with the outbursts 
in the two systems is given in Table 1 . The behavior of the 
sequences of outburst durations t& and quiescent intervals t q 
during a supercycle are quite different. The tb sequence shows 
an increase between SOs in both systems. The t q values re- 
main relatively constant between SOs in V1504 Cyg, whereas 
in V344 Lyr they increase to a local maximum roughly half 
way between SOs before returning to their initial small value 
by the time the next SO begins. Thus the t q values are roughly 
continuous across SOs. In both systems the variation of t q 
values between SOs appears to be inconsistent with the pre- 
diction of the thermal-tidal instability model, which shows a 
monotonic increase in t q between SOs (Ichikawa, Hirose, & 
Osaki 1993, see their Figure 1; Osaki 2005, see his Figure 3). 

The asymmetry of the NOs is similar in the two systems, 
f r ~ 0.2 - 0.3 in VI 504 Cyg and f r ~ 0.25 - 0.35 in V344 
Lyr. According to theory, the degree of outburst asymmetry 
reveals whether the outburst is triggered at small or large radii 
in the disk (Smak 1984; Cannizzo, Wheeler, & Polidan 1986. 
hereafter CWP; Cannizzo 1998, see his Figure 6 and asso- 
ciated discussion). Outbursts triggered at large radii tend to 
have fast, asymmetric rise times (/ r < 0.5), while those trig- 
gered at small radii are more symmetric (/ r ~ 0.5). Smak 
referred to these two types as Type A and Type B. respec- 
tively, whereas CWP used the terminology “outside-in" and 
“inside-out" outbursts. The low f r values in both systems 
studied in this work are consistent with the Type A/outside-in 
outbursts. In is interesting that the one outburst which seems 
to be certifiably symmetric and therefore presumably inside- 
out — outburst 1 of V344 Lyr — occurred after an unusually 


long quiescence interval, which may have led to a more filled 
disk at the end of quiescence (compare this outburst with that 
shown in Figure 5b of Smak 1984 and Figure 3 of CWP). 
During this part of the V344 Lyr light curve there were strong 
negative superhumps (Wood et al. 201 1 , see their Figure 1 3), 
indicating a tilted accretion. This would have the effect of 
spreading the accretion stream impact over a large range in 
radii, thereby lengthening the time to build up to the next out- 
burst. 

The mean durations of SOs for VI 504 Cyg and V344 Lyr 
are 11.9 d and 16.8 d, and the e— folding rates of the (slow) 
decay during SO are 7.9 d mag” 1 and 12.0 mag -1 , respec- 
tively. Thus the SOs span 1.5 e— foldings in VI 504 Cyg, and 
1.4 e— foldings in V344 Lyr. Given that the SO accretion disks 
are close to steady state, the depletion fraction f d in the stored 
disk mass during SO for the two systems, 

f d ~ 1 - e” 4(SO)/td(so) , (1) 

is ~0.78 and ~0.75, respectively. (These numbers probably 
over-estimate the depleted mass fraction because they encom- 
pass portions of the SO fainter than the e— folding plateaus.) 
These depletion fractions are larger than the fd — 0.2 inferred 
for long outbursts in DNe above the period gap (e.g., Can- 
nizzo 1993), and stem from the larger a bot /a co id ratios in the 
SU UMa systems. The ratio t d (SO)/td(NO) ~ 20 is consis- 
tent with the previous interpretation of fast decays being due 
to the action of a cooling front within the disk, and slow de- 
cays being the result of gradual accretion onto the WD, and 
therefore governed by the viscous time scale in the outer disk. 
The viscous and thermal time scales can be written in terms 
of the local Keplerian rotation frequency fi and the local disk 
height to radius ratio h/r (Pringle 1981). The viscous time 
scale 

1 /r\ 2 

u = (ft) ( } 

and the thermal time scale 

tth = a hot Q (ft) (3) 

differ by a factor h/r, therefore h/r « 1/20. 

C 10 present an equation for the rate of decay of a SO, taken 
from Warner (1995a, 1995b), which is based on the viscous 
time in the outer accretion disk 

t d « t„ = 17 d Ph /4 m 1 ^, (4) 

where cthot,-i = ahot/0.1 and mi = M\/Mq. This 
e— folding time scale should approximate the decay rate, ex- 
pressed in d mag” 1 . Hence in a large sample of DNe the 
prediction is that the rate of decay during SO should vary 
as P/” 5 . Our analysis shows, however, that in going from 
V1504 Cyg with P h = 1.67 to V344 Lyr with P h = 2.1. the 
SO decay rate increases from 7.9 d to 12 d. One does expect 
object-to-object variations in NO and SO properties which 
can be large compared to overall, secular trends with orbital 
period, but given the dependencies in Equation (4) and the in- 
ference from a large number of DNe that cthot — 0.1, there is 
no way to account for such a steep variation, t d 3c P/' 8 . 

It is worth pausing for a moment and looking at how the 
SO slow decay rates for V 1504 Cyg and V344 Lyr fit into the 
larger picture. Table 2 lists the slow decay rate of SOs for SU 
UMa stars taken from the literature, most of which are derived 
from VSNET data (see Kato et al. 2009, 2010). We only show' 
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TABLE 2 

Slow Decay of Superoutbursts in SU UMa stars 


System 

decay rate (d mag l ) 

Ph 

Reference 

SW UMa 

10 . 

1.36 

1 

WX Cet 

10 . 

1.40 

2 

FLTra 

7.7 

1.41 

3 

VI 028 Cyg 

8.3 

1.45 

4 

HO Del 

7.1 

1.50 

5 

GO Com 

6.7 

1.58 

6 

NSV 4838 

10. 

1.63 

7 

VW CrB 

10.3 

1.70 

8 

EG Aqr 

7.1 

1.83 

9 

V503 Cyg 

9.1 

1.86 

10 


References. — 1. Soejima et al. (2009), 2. Kato et 
al. (2001), 3. Imada et al. (2008b), 4. Baba et al. (2000), 
5. Kato et al. (2003), 6. Imada et al. (2005), 7. Imada et al. 
(2009), 8. Nogami et al. (2004), 9. Imada et al. (2008a), 10. 
Corrado et al. (2000) 


data for the VW Hyi subclass for which the slow decay rate is 
relatively constant over the course of the SO. In Figure 1 1 we 
plot the SO decay rates for the systems shown in Table 2, as 
well as V1504 Cyg, V344 Lyr, WZ Sge and U Gem. For WZ 
Sge we only plot the value for the initial, steep rate of decay, 
4 d mag -1 . As noted in the Introduction, the SO decays in 
the WZ Sge stars exhibit a strong deviation from exponential- 
ly presumably due to the accretion of virtually all the stored 
mass. As the disk mass decreases, one expects a concomitant 
increase in the viscous time, leading to a pronounced concave 
upward SO decay. Therefore only the very initial SO decay 
rate gives a valid representation for the viscous time within 
the disk corresponding to the stored mass. As regards DNe 
above the period gap, virtually none of the long outbursts ever 
observed exhibited enough dynamic range in their overall de- 
cay that one could reliably extract the rate of decay. Only for 
the unusually long, 1985 October outburst of U Gem could 
this determination be made (Cannizzo, Gehrels, and Mattei 
2002), hence the single data point above the period gap in 
Figure 1 1 . If one takes the overall t d — P/, trend between WZ 
Sge and U Gem as a true physical baseline, t d oc P/ 6 , the 
trend between V1504 Cyg and V344 Lyr in logt d — logP^ 
space, t d oc P/ 8 , is roughly the same, although the normal- 
ization is offset slightly. Given the vertical scatter in Figure 
1 1 for the other SU UMa systems, this coincidence may be 
fortuitous. 

The failure of Equation (4) to account for the difference 
between the SO slow decay rate in V1504 Cyg and V344 
Lyr indicates that a more sophisticated expression may be re- 
quired. One obvious deficiency is that there is no dependence 
on a co id which controls the available amount of fuel in qui- 
escence to power a SO. The physics behind the viscosity in 
quiescence is unknown at present, and one could in principle 
have large variations in a C oid between systems at similar or- 
bital periods. Let us go back to the definition of the viscous 
time scale t v = (ahot.^)~ l ( r A) 2 - Since the controlling time 
scale in this situation is the slowest one. this expression is to 
be evaluated at the outer disk radius r d . Following Cannizzo 
(2001b), we consider the hot state of the disk bv adopting a 
Rosseland mean opacity k — 2.8 x 10 24 g cm and 
mean molecular weight p = 0.67. We can then write the vis- 
cous time as 


>—r i oj 14 13/14 
5-7 dm/ r' c 


—8/7 
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Fio. 6.— The values of the outburst durations fj, for the two systems, taking 
a cut line Kp cat = 16 for V1504 Cyg and Kp cu t = 17.9 for V344 Lyr. 
The vertical dotted lines ( red) indicate the superoutbursts. The values of tj, 
for the superoutbursts are t;,(d) = 12.0, 1 1.9, 1 1.8, and 11.9 for V1504 Cyg, 
and f;,(d) = 16.6, 17.1, 16.6, and 16.9 for V344 Lyr. 

where m x = M WD /lM e , r d> i 0 = r d /( 10 10 cm), and E (r d ) 

is the surface density at the outer disk edge. 

Looking at this expression, we see that the Warner ( 1 995a, 
1995b) equation already contained an implicit assumption 
about E, as it only contains cthot and P d . To derive a more 
meaningful expression for t d « t v in the context of SOs we 
need to take into account two additional pieces of information, 
namely that during a SO essentially all the mass accumulated 
during quiescence is accreted onto the WD, and thereby lost 
from the disk, and also that the accumulated mass depends on 
the critical surface density E max , which depends in turn on 
ctcoid ■ This scaling is given by (Lasota 2001) 

^max — io.4 g Cm Tjq . (6) 

Since the attainment of E max is the triggering mech- 
anism for a DN outburst, the quiescence disk mass 
~(7r/2)r^E max (r < i, a C oid) at the time of SO onset can be 
written as 

AM cold » 4.2 x 10 21 g mr 0 - 38 r^ 4 Q-°f . (7) 

Once the SO has started, the mass that was present in the 
non-steady state quiescent disk gets redistributed into a quasi- 
steady disk with M (r) ~ constant, and therefore with a radial 
profile E(r) oc r~ 3 '' 4 (SS). Therefore the mass of the hot disk 
can be written as 


AA/ h 


27 rre(rE(r, a hot ). 


(8) 


We now set E (r, c*hot) equal to its standard functional form 
for a SS disk with gas pressure and free-free opacity, Ehot = 

E 0 a/// J r 1 ” 0 3/,4 . This is the same form that entered into Equa- 
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Fig, 7. — The values of the asymmetry parameter f r = t r i se /4 for the 
normal outbursts in the two systems. Only values for which the times of 
outburst maxima can be reliably determined are shown. The same cut line 
values are used as in Figure 6. 


tion (5). Equating the disk mass at the end of the quiescent 
state with that at the beginning of the SO, i.e„ AM co id — 
AM hot , gives 

5 

^hot ir d) “ Yg ^max (U), (9) 

the relevant E value for Equation (5). Therefore we can 
rewrite Equation (5) as 

U = 80.2 d m? ' 52 etcofd OfTot r d ; to ■ 00) 

This expression differs from Equation (5) in that it contains 
ct co id, reflecting the fact that the total amount of matter stored 
in quiescence enters into the viscous time scale in the resultant 
SO disk. 

In extreme mass ratio systems such as the SU UMa stars, 
the Roche lobe of the primary asymptotes to ~0.82 of the 
binary separation (Eggleton 1983). Also, the relation be- 
tween orbital separation and orbital period is aio = 3.53 P h , 
where aio is the orbital separation in units of 10 10 cm. If we 
further take the outer disk radius to lie at 75% of the pri- 
mary Roche lobe, we have r<uo = 2.19 Pf/ 3 - Since it is fairly 
certain that a in outburst is always close to 0.1, a+ot been 
scaled to that value in the previous equations. On the other 
hand, the a value in quiescence is uncertain, and has been tra- 
ditionally adjusted by trial and error for each system. For SU 
UMa stars, workers have typically had to adopt much smaller 
a co id values than for DNe above the period gap. Given this 
uncertainty, let us take a co id to be a power law in orbital pe- 
riod. a C otd = a c ,o Pjf • Making these substitutions into Equa- 



FlG. 8. — The values of the quiescence durations t q between consecutive 
outbursts for the two systems. As with Figure 6, the vertical dotted lines (red) 
indicate the positions of the superoutbursts. The same cut line values are used 
as in Figure 6. 

tion (10) we obtain 


t u = 113 d to?' 52 ' 


a, 


0.36 

c,0 


a 


-1.14 
hot,— 1 


p 0. 29+0. 36^ 

J h 


(ID 


The dynamic variation associated with one magnitude (2.512) 
is 0.924 times that associated with one e— folding (2.718). 
Therefore if we associate the rate of SO decay, expressed in 
d mag -1 , with the viscous time at the outer edge, which is an 
e— folding time, then Equation (11) needs to be corrected by 
a factor 0.924. 

By comparing the SO decay rates of V1504 Cyg and V344 
Lyr, we infer a decay rate fd = 0.924£„ = 12 d ( P ^/ 2 . 1) 1 ' 8 = 
3.16 d Pi' 8 . Furthermore, adopting rn x = 0.6 for the SU 
UMa systems gives to ?-° 2 = 0.767, and therefore 

80.2 d a 0 c 36p0.29+0.36i> = 3 d pl.S ( 12 ) 


Solving for the constant and exponent in the a co id scaling 
gives a c , 0 = 1.25 x 10 -4 and ip = 4.2 — a steep depen- 
dence of a co id on orbital period. For V344 Lyr at Ph = 2.1 
we obtain a cold ~ 2.8 x 10 -3 , which accounts for why CIO, 
in their modeling of V344 Lyr, had to adopt the small value 
Q+oid = 0.0025. The inference ip = 4.2 is only valid, how- 
ever, if the trend between V1504 Cyg and V344 Lyr were 
representative of a trend with P/,. which seems unlikely given 
the data for the other SU UMa systems. Based only on the SU 
UMa data, it would be more likely that a c ,o has a large vari- 
ation between the two systems. As noted earlier however, if 
one takes the WZ Sge <-* U Gem interpolation to be physically 
meaningful, then the V1504 Cyg «-» V344 Lyr interpolation 
appears to be fortuitously coincident. 

A steep dependence of a- co id on P/, may have observational 
support from another line of argument. In a recent study of the 
aggregate properties of DNe. Patterson (201 1. see his Figure 
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FIG. 9. — A composite of four superoutbursts in V1504 Cyg, positioned in 
time so that t = 4 d corresponds to the crossing of Kp = 16 during super- 
outburst onset. The line segments ( red) indicate the linear least squares fit, 
for the range of data indicated. The rate of decay for the four superoutbursts 
are as shown in each panel. The mean value is 7.9 d mag - 1 . 

1 1) finds that the recurrence time for SOs has a steep inverse 
relation with binary mass ratio f r ecur(SO) oc g-2.63io.i7, 
where q = M2 /Mi- Since the secondaries in cataclysmic 
variables are, by definition, in contact with their Roche lobes, 
this implies that M 2 scales with orbital period, and hence 
frecur(SO) oc p- 2 ' 63±0 - 17 . Given that (1) SOs accrete vir- 
tually all their stored mass AiW co id, (2) the mass feeding rate 
into the outer disk from the secondary Mp is thought to be 
relatively constant at ~10 -10 Mq yr _I for systems below the 
period gap (Kolb et al. 1998; Howell et al. 2001), and (3) the 
recurrence time (Cannizzo et al. 1988) 


^A^cold ~ 3.14 - 

lilfrl X,d c 


it is quite surprising that PecurfSO) apparently has a strong 
inverse relation with orbital period. As one progresses to 
shorter orbital period, one expects more frequent SOs since 
the disks are smaller and the critical mass is reached more 
quickly. We are seemingly drawn to the same conclusion as 
previously, that a co id must have a strong dependence on or- 
bital period. Howell et al. (2001; see their Figure 2, panel (a)) 
indicates a rough dependence dlog M;- /dlog l\ a 1 for 
the SU UMa stars. If we again set a co id oc P / , substitute Ph 
for r ( i in Equation (13), and solve for ip, we obtain ip = 4.5. 

If the overall trend in t d for long outbursts between WZ Sge 
and U Gem, which is roughly matched (coincidentally) by that 
between VI 504 Cyg and V344 Lyr, is real, then two indepen- 
dent lines of reasoning imply a very steep dependence of a co id 
on orbital period. Why should one have Qhot — constant. 


FlG. 10. — A composite of four superoutbursts in V344 Lyr, positioned in 
time so that t = 3 d corresponds to the crossing of Kp = 18 during super- 
outburst onset. The line segments ( red) indicate the linear least squares fit, 
for the range of data indicated. The rate of decay for the four superoutbursts 
are as shown in each panel. The mean value is 12 d mag -1 . 


40 


r- 20 

1 

hfl 

cd 

10 
8 

3 6 

T3 

^ 4 


1 2 3 4 5 

P 
r h 

Fig. 1 1 . — A comparison of superoutburst slow decay rates of other SU 
UMa stars, from Table 2, with VI 504 Cyg and V344 Lyr. Also shown are 
WZ Sge (initial superoutburst decay rate only) and U Gem. The vertical 
dashed lines indicate the 2 — 3 hr orbital period gap for DNe. 

while at the same time Q co id oc P d , where r ~ 4.2 - 4.5? 
Current thinking on the physical mechanism for viscosity in 
accretion disks centers on the magnetorotational instability 
(MRI - Balbus & Hawley 1998), an instability based on the 
shearing amplification of a weak magnetic field. The non- 







linear saturation limit of the MRI in ionized gas is thought 
to explain why «hot — >~0.1 (McKinney & Narayan 2007a, 
2007b). However, as the disk cools to the quiescent state and 
the supply of free electrons drops precipitously due to the ex- 
ponential dependence of partial ionization on temperature be- 
low ~10 4 K, the magnetic field no longer couples effectively 
to the gas, and fails to provide an effective viscosity. Menou 
(2000) tried to model the complete accretion disk limit cycle 
model for DNe using only the MRI, and failed. He discussed 
the failure of the MRI in quiescence, and predicted that if tidal 
torques (a completely different physical mechanism from the 
MRI) dominate the angular momentum transport in the qui- 
escent accretion disks in the SU UMa stars, one would ex- 
pect an anticorrelation between f r ecur(SO) and q. He presents 
data for six SU UMa stars showing an anticorrelation. (Pat- 
terson [2011] shows the anticorrelation in much greater de- 
tail.) Menou notes that theories of tidally induced spiral waves 
or shock waves indicate a strong correlation between effec- 
tive torque and q (Papaloizou & Pringle 1977; Goldreich & 
Tremaine 1980). Our finding a co id oc P^' 2 seems to support 
the view that tidal torquing is the dominant angular momen- 
tum transport mechanism operating in quiescent accretion 
disks in interacting binaries. 

4. CONCLUSION 

We have analyzed long term Kepler light curves of VI 504 
Cyg, containing 40 outbursts, and V344 Lyr, containing 38 
outbursts. Each system showed four superoutbursts. The find- 
ings are: 

(1) The decays of the NO are faster-than-exponential, with 
most of the deviation occurring near maximum light. Near 
quiescence the decays are close to exponential, ~0.7 d mag" 1 
for VI 504 Cyg and ~0.6 d mag -1 for V344 Cyg. These rates 
are in line with that expected from the Bailey relation. Our 
values are based on frequency histogram distributions of only 
~30— 35 NOs and should improve as more outbursts are ac- 
cumulated. 

(2) The outburst durations % for the NOs increase by a fac- 
tor ~ 1.2— 1.7 between consecutive superoutbursts. 

(3) The quiescent intervals t q between normal outbursts 
remain relatively constant for three of the four supercycles 
in VI 504 Cyg, whereas in V344 Lyr the t q values increase 
to a local maximum about half way between superoutbursts. 
There was also an anomalously long quiescence interval at 
the start of the V344 Lyr light curve, potentially the result of 
a tilted disk. Neither of the observed t q trends are consistent 
with Osaki's thermal tidal instability model for SOs, in which 
one expects a monotonic increase in t q values accompanying 
the monotonic build-up in disk mass and angular momentum 
leading to a SO. 

(4) The asymmetric shape of the NOs is coincident with the 


triggering for the outbursts being outside-in. 

(5) The inference of the steep relation a co id oc Pjl gives 
strength to the notion of tidal torquing as the dominant an- 
gular momentum transport mechanism operating in quiescent 
accretion disks in interacting binaries. 

If one considers the slow decay rate of the long, viscous 
outbursts spanning DNe from WZ Sge to U Gem, one in- 
fers an overall variation oc -P/J' 6 - Given the scatter in the 
decay rate for the 10 well-studied SU UMa systems shown 
in Table 2, it may be fortuitous, given the general scatter in 
DN outburst properties at a given orbital period, that the de- 
cay rate between VI 504 Cyg and V344 Lyr shows a similar 
law oc P^' 8 . These dependencies are much steeper than the 
P °' 25 expected from the previous theoretical expression, due 
to Warner (1995a, 1995b). By starting with the definition for 
the viscous time scale in the outer disk during superoutburst, 
and taking into account the scaling for material accumulated 
during quiescence, we derive a more physically motivated for- 
mula, and find that a steep dependence of the quiescent state 
value of a is required, a co id oc P^' 2 . In simple terms, the dra- 
matic increase in stored quiescent mass with decreasing or- 
bital period brings about a decrease in the viscous time scale 
of matter in the superoutburst disk, which is formed out of re- 
distributed gas from the quiescent state. Recent work by Pat- 
terson (201 1) on the recurrence time scale in SU UMa stars 
for SOs also supports a steep dependence of a co id on P/,. 
The long outbursts in systems above the period gap and su- 
peroutbursts in systems below the gap both represent viscous 
decays, with the absence of transition waves in the disk. It is 
somewhat counterintuitive that the recurrence time for these 
long, viscous outbursts varies as P ; " 2,6 , and yet their decay 
time varies as P^ 8 . The steepness of the a co id (Ph) relation 
explains why accretion disk modelers have had to adopt such 
small a co M values for the SU UMa stars, compared to DNe 
above the period gap, and also probably accounts for the fact 
that the extreme mass ratio black hole transient systems un- 
dergoing limit cycle accretion disk outbursts, like A0620-00, 
have such long recurrence times for outbursts. 

The Kepler data have already provided important con- 
straints on the physics of accretion disks, in particular the 
Shakura-Sunyaev a parameter: CIO showed that an interpo- 
lation between a co id and exhot based on the degree of partial 
ionization of gas was needed to account for the shoulders in 
the SOs of V344 Lyr, and in this work we have shown that the 
scaling of SO decay rate with orbital period mandates a steep 
dependence of cx co id on orbital period, enforcing the notion of 
tidal torquing as being the dominant viscosity mechanism in 
quiescent accretion disks. 

We acknowledge the contributions of the entire Kepler 
team. 
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